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Abstract 



We derive the Ginzburg-Landau free energy functional of color su- 
perconductivity in terms of the thermal diagrams of QCD in its pertur- 
bative region. The zero mode of the quadratic term coefficient yields 
the same transition temperature, including the pre-exponential factor, 
as the one obtained previously from the Fredholm determinant of the 
two quark scattering amplitude. All coefficients of the free energy can 
be made identical to those of a BCS model by setting the Fermi velocity 
of the latter equal to the speed of light. We also calculate the induced 
symmetric color condensate near and find that it scales as the cubic 
power of the dominant antisymmetric color component. We show that in 
the presence of an inhomogeneity and a nonzero gauge potential, while 
the color-flavor locked condensate dominates in the bulk, the unlocked 
condensate, the octet, emerges as a result of a simultaneous color-flavor 
rotation in the core region of a vortex filament or at the junction of 
super and normal phases. 
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1. Introduction. 
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The properties of hadronic matter under unusual conditions, high temperature or 
high baryon density, have been the subject of intense theoretical investigation, especially 
during the last few years. Relativistic heavy-ion collision experiments probe the quark- 
gluon plasma both at nonzero temperature and nonzero chemical potential. Quark matter 
at high baryon density and low temperature is expected to exist inside neutron stars. 
The exploration of these exotic states of matter will hopefully shed light on many non- 
perturbative aspects of quantum chromodynamics. 

Although great progress has been achieved in our understanding of QCD along the 
temperature axis in the two-dimensional phase diagram of QCD with respect to temper- 
ature and chemical potential, based primarily in our ability to perform simulations on 
the lattice, the problem of low temperature, high density hadronic matter is less well un- 
derstood mainly because the fermion sign problem at nonzero chemical potential makes 
numerical simulations practically impossible at high baryon density. Thus, many issues of 
the physics near the chemical potential axis can only be addressed analytically. 

At ultra high chemical potential, n » AqcD) the interaction between two quarks 
near their Fermi surface is dominated by one gluon exchange because of the asymptotic 
freedom. Such an interaction can be decomposed into a color symmetric channel and a 
color antisymmetric channel. For a pair of quarks within the color antisymmetric channel 
propagating in opposite directions, both color-electric and color-magnetic parts of the 
interaction are attractive. Therefore the Fermi surface is unstable against the formation of 
Cooper pairs and the quark matter becomes superconducting below a certain temperature 
[1]. Such a color superconducting state is expected to survive at medium high chemical 
potential, fj, > Aqqd (a* = 400 ~ SOOMeV), which corresponds to the baryon density 
achieved in RHIC and in the core of a typical neutron star. The instanton interaction may 
also contribute to the condensation at medium high chemical potential [2]. 

Theoretically, there are two approaches which explore the physics of the color super- 
conductivity, one that utilizes an effective action involving only quark degrees of freedom 
while the other uses the fundamental QCD action. Using the first approach, important 
progress has been made by Alford, Rajagopal and Wilczek, who recognized that the color- 
flavor locked condensate represents the minimum of the bulk free energy and thus charac- 
terizes the super phase [3] . As this condensate possesses all the symmetries of low energy 
hadronic matter, the cross-over from one phase to the other is likely to be continuous [4]. 
The second approach [5], [6], [7], [8], [9], [10], [11], [12] yielded as a main result an asymptotic 
formula of the superconducting energy scale in the weak coupling limit. For example, the 
transition temperature with Nc colors and Nf flavors is given by 
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with g being the QCD running couphng constant at /u, and 7 = 0.5772.. the Euler constant. 
The non-BCS behavior of the exponent stems from the long range propagation of the mag- 
netic gluon being exchanged and the 0(1) term of the exponent is due to the suppression 
of the quasi-particle weight at the Fermi level, another consequence of the long range mag- 
netic interaction. Other contributions which are subleading in g [11], including the scheme 
dependence of the ultraviolet renormalization do not appear in (1.1). Remarkably, the 
formula (1.1), even when is extrapolated to the medium high chemical potential, gives rise 
to a very small ratio ksTc/ jJ-- 



In this paper, we shall derive the Ginzburg-Landau free energy functional of color 
superconductivity at ultra high chemical potential. The conventional method of derivation 
was developed by Gorkov [13], who generalized the superconducting gap equation to a set of 
equations for ordinary and anomalous Green functions, the latter being appropriate to an 
inhomogeneous system. Although this formalism is easy to implement in the case of a four- 
fermion interaction it becomes rather cumbersome when one uses the fundamental QCD 
action with an energy dependent gap. Furthermore the systematics is not transparent. 
Instead, we shall employ the formulation developed by Jona-Lasinio [14], in which the 
Ginzburg-Landau free energy is identified as the generating functional of proper vertices. 
Each coefficient of the Ginzburg-Landau free energy is associated with a set of thermal 
diagrams in the normal phase. As we shall see the transition temperature corresponds to 
the zero mode of the inverse di-quark propagator and this criterion is rigorously equivalent 
to the one used in [10], [11], in which the transition temperature was identified with the 
zero of the Fredholm determinant of the di-quark scattering kernel. For a homogeneous 
system, the derived Ginzburg-Landau free energy agrees with the recent result of lida and 
Baym [15] to the leading order in g and in the absence of the color symmetric component 
of the condensate. The color-flavor locked condensate remains the energetically favored 
minimum in the bulk near Tc- Including the inhomogeneity, the Ginzburg-Landau free 
energy for QCD agrees exactly with that of a BCS superconductor upon setting the Fermi 
velocity of the latter to the speed of light, despite the non-BCS scaling of the transition 
temperature (1.1). 



Although the one gluon exchange is repulsive in the color symmetric channel, a small 
component of the condensate in this channel was found to be induced at T = [3], [16], 
[17]. We shall analyze this mechanism as T approaches Tc from below. We find that the 
magnitude of the induced color symmetric condensate scales like the cube of that of the 
color antisymmetric condensate. The contribution to the free energy at equilibrium from 
the former is suppressed relative to that from the latter by a factor of the order of ^^^"'^^ , 
which is larger than the one that was speculated in [18]. 



As an important application of the Ginzburg-Landau theory, we shall also examine 
the coupling between the color-fiavor locked condensate and the unlocked one through a 
gauge potential for an inhomogeneous system. We find the inevitable appearance of the 
unlocked condensate in the core region of a vortex filament or at the border of the super 
phase while the locked condensate remains dominant in the bulk. 
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This paper is organized as follows: In section 2 we review the Jona-Lasinio method 
which we shall employ in this paper while in section 3 we discuss different characheristics 
of the order parameter of color superconductivity, which corresponds to a non-zero vacuum 
expectation value of a quark bilinear. Next in section 4 we calculate the Ginzburg-Landau 
free energy functional to the quadratic order in the order parameter and in section 5 we 
extend the calculation to the quartic order. Finally in section 6 we introduce a non-zero 
total momentum to the diquark condensate and the corresponding gauge coupling and 
examine the application to the realistic case with three colors and three flavors. 

2. General Formulation. 



In this section, we shall review the Jona-Lasinio method [14], which we shall em- 
ploy in the microscopic derivation of the Ginzburg-Landau free energy functional of color 
superconductivity through out this paper. 

Let £s be the Lagrangian density which describes a field theoretic system in Euclidean 
space. The free energy of the system at temperature T, F{T) can be expressed in terms 
of a functional integral 



where /? = l/ksT, Se{(I>) = Jq dr J d^rCE{(j>) and 4'a stands for the collection of the 
field variables of the system with the subscript a denoting its spacetime and internal 
coordinates. Periodic boundary conditions (anti-periodic) in Euclidean time r are imposed 
on the bosonic fields (fermionic). 

The order parameter of the system corresponds to the nonzero expectation value of 
an operator, which may be one of the elementary fields of the theory or a composite 
operator and realize a nontrivial representation of the symmetry group of the system. To 
explore the stability of the system against the order parameter, we first trigger it with an 
external source, i.e. 




a 



(2.1) 



Se{4>) ^ Si{4,) = SE{ct>) + 5]( J:$a + JaK)- 



(2.2) 



a 



The corresponding shift of the free energy, W{J), is given by 



e 



(2.3) 
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and can be expanded in powers of J and J*, i.e. 

^(•^) = -lT.^abJ:Jb + ^ E GabcdKHJcJd + 0(1 J|«), (2.4) 
ab abed 

where each coefficient corresponds to a set of connected diagrams that contains a fixed 
number of J's. In particular, T>ab corresponds to the propagator of The order param- 
eter triggered by J reads 

Ba =< $a >-(^JJZ = - J2'^'^bJb + J2'^'^bcdJ;'^cJd + 0{\J\^). (2.5) 



After performing a Legendre trasformation, 



r{B) = W{J) - Y^iJlBa + JaBl), (2.6) 

a 

we find that 

^ -Ja. (2.7) 



SB, 

Inverting the series (2.5), we obtain that 



r{B) = J2 MabBlBk + \Y. QabcdBlBlB.Bd + 0{\B\'') (2.8) 



2 

ab abed 



with 
and 



Mab = {T>-%b (2.9) 

Gabed= Yl ■^la'J^bb'Ga'b'c'd'Mce'Mdd'- (2.10) 



a'b'c'd' 



The functional r(S) is nothing but the Ginzburg-Landau free energy functional of the order 
parameter. When the external source is removed, the minimization of F at a nonzero 
B signifies the condensate phase. For a second order phase transition, the sign of the 
quadratic term plays a decisive role. 



Since the quartic term is positive and we are interested in the physics near Tc we shall 
subsequently neglect higher terms. Few remarks are in order: 



1. If the operator is an elementary field or a quadratic function of the elementary 
fields, the manipulation leading from (2.2) to (2.8) can be carried out for a free field system 
with all the coefficients of (2.8) known explicitly. In this case F is always positive for a 
nonzero B and the only minimum of F occurs at S = 0. This is also the case with a 
repulsive interaction. 
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2. The presence of an attractive interaction makes a difference. Let us diagonalize the 
matrix A4 and expand the order parameter B in terms of the cigcnmodes of Ai. As we shaU 
see below, in the presence of a Fermi sea, no matter how weak is the attraction, there is 
always a mode, whose corresponding eigenvalue becomes negative below a sufficiently low 
temperature Tc- We shall refer to this mode as the pairing mode in contrast to the others 
(whose eigenvalues are positive) which we shall refer to as non-pairing modes. Technically, 
we may switch off the interaction when dealing with the non-pairing modes in the leading 
order weak coupling approximation since the interaction is merely perturbative to them. 
The same approximation applies to the quartic term of (2.8), since the pairing singularity 
of Gabcd is cancelled by the zero of Al's upon amputation. 



3. It is not always true that the order parameter which minimizes F consists of only 
the pairing mode. In the presence of a term that is cubic in the pairing mode and linear 
in a nonpairing mode, an expectation value of the nonpairing mode will be induced upon 
minimization, the magnitude of which will be proportional to the cubic power of that of the 
pairing mode. As the critical temperature Tc is approached from below, the pairing mode 
condensate scales as {Tc — T) 2 and the induced nonpairing mode scales as {Tc — T) 2 . As 
we shall see below, this is the case with the sextet condensate in the color-flavor locked 
phase of three colors and three flavors. 



4. As is the case for all mean field theories in three dimensions, the weak coupling 
calculation can not be extended to the temperature arbitrarily close to the critical tem- 
perature Tc, since the perturbative expansion in the coupling constant is not uniform in 
T — Tc- The size of the region where our approximation fails, the critical region, can be 
estimated according to the Ginzburg criterion. Because of the agreement of our derived 
Ginzburg-Landau free energy with that of a BCS model, we conclude that the width of 
the critical window is given by 



IT -Tel fkuTc 



Tc \ n 



(2.11) 



which is exceedingly narrow. Technically, the critical point prevents us from continuing 
r{B) obtained in the normal phase to the super phase at J = 0. A nonzero triggering 
source J is required for a smooth cross over to the super phase. To be more specific, 
the normal phase F functional can be carried over to the super phase following the path 
ABCD in Fig.l. 



In the appendix A, we shall utilize the Jona-Lasinio formulation in order to derive 
the Ginzburg-Landau free energy functional for a BCS model as an illurstration and as 
comparison to the analogous expression for color superconductivity. 
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3. The Order Parameters of QCD at High Baryon Density. 



Consider now an SU (Nc) color gauge theory coupled to Nf flavors of massless quarks 
at temperature T and chemical potential The Eucliden action of the system reads 



= J^r J dV[iF^,F^,-V'7M(^-^^^M)V' + MV'V + ---], (3.1) 



where Fj^^ = a^Al, - a^^Aj^ +5f/^"^A^A^, = A^^T'' with being the SU{Nc) generators 
in the fundamental representation and the dots represent gauge fixing and ghost terms. 
Since the action (3.1) is diagonal with respect to both flavor and chirality, the corresponding 
indices have been dropped in (3.1). Our gamma matrices are all hermitian with 75 = 
71727374- 

The fermion field can be decomposed into eigenfunctions of 75, which in turn can be 
separated into the positive energy and negative energy parts 

rfir.r) = rfL+{r,T) + rfL-{r,T) + rfR+{r,r) + rfR-{r,r), (3.2) 

where we have restored the color-flavor indices for clarity. The Fourier expansion of 
^L{R)±i^^^) is given by 





1 


p 




rfL-ir,r) = - 


1 


p 




rfR+{r.T) = - 


1 


p 


-'^-^alp^^URi^ 


^'fL-{f,r) = - 


1 


p 


-'"-^blp^RM-P) 



(3.3) 



and parallel expressions in terms of p ^, b'j p j^^, cifpji and Iffpji exist for '^''f L{R)±i'^i ''')' 
where P = (p, — fn) denotes an Euclidean four momentum with z/„ = 2mTkBT (n = 
±|, ±|, ...) the Matsubara energy. The four component spinors UL(fi-j{p) satisfy the Dirac 

equation {'^aP — '^1 ■ 'P)ul{r){'P) = and are normalized according to ^t^^.^^|(^?)•ux,(fl;)(^?) ~ 1. 
Furthermore, we have '^5ULi'P) = ul{p} and '^5Ur{'p) = —ur{j>). Similarly fL(i?)(p) obeys 
the Dirac equation and the chirality conditions. The symbol stands for the volume of 
the system and as f2 1— > 00 
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To trigger the long range order, we add to (3.1) a source term, 

^S=l Yl f dr, I d'n f dT2 I d'r, fe,,,,(fi, n; f2, r^) 

h=L,R;s=±'^^ J Jo J ^ 

(3.5) 

with 

'^hh,h,siri, -n; f2, -T2) = j;X.,,(ri, n; r2, r^) (3.6) 

and 

J}'j:,HAri,ri;r2,r2) = J-;,V^,V,(r-2, r^; n, n). (3.7) 

The condition (3.6) is the Wick rotation of the complex conjugation with a real time and 
the symmetry property (3.7) is evident from the anticommutation relation of the ip^s and 
the antisymmetry of a2- Both of them are shared by the induced order parameter 

mihA^^^^^-^^^^^^) ^< rA,,Jr,,n)a2r/,,,,sir2,r2) > . (3.8) 

To trigger a homogeneous condensate, the source should depend only on the relative Eu- 
clidean coordinates, (ri — f2,Ti — T2) and its Fourier expansion reads 

JC1C2 f \ _ Jp-{ri-r2)-il^ri{Tl-T2) JC1C2 (p\ /o q^ 

*ith J%t,h,s{P) = JZt,h,si-P)- Correspondingly, 

with -Byjjj h si-P) ~ PfT/i h s(~P')- '^^^ triggering part of the action (3.5) becomes 

= T!iJ^%,>^Ana%,.paf^,p + jrji,HAnd%,pd%,.p] (3.11) 
p 

with the sum over P extending only to the half p^space and J = J*. 

As the condensate consisting of quarks of left helicity contributes identically with the 
condensate of quarks of right helicity the total Ginzburg-Landau free energy functioncal 
of the system can be written 

TtotaiiBL, Br) = T{Bl) + TiBn) + ■■■ (3.12) 

where the • • • indicate terms which include both Bl and Br and arise due to instanton 
effects [19], [20]. Such terms are exponentially small in comparison with the first two terms 
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of (3.12) but nevertheless contribute to fix the relative phase between Bl and Br. Since 
the even parity is favored for the ground state, we expect that 

^n,LA^i^rur2,r2) = B}l%^j,^^i^^^^ (3.13) 

and (3.12) becomes 

Ttotai{BL,BR)^2T{B) (3.14) 

where we have omitted the extra terms in (3.12). T{B) is the Ginzburg-Landau free energy- 
functional with a single helicity. It's quartic expansion 

T{B)=T2{B)+T4{B) (3.15) 

will be evaluated in the subsequent sections with the helicity subscripts L, R omitted for 
clarity of notation. 

At high baryon density, the pairing correlation extends hardly to the states below the 
Dirac sea and we may set J71 = with J7+ = JT". 

In this paper, we shall consider only the s-wave pairing for which both Jflfl (P) and 
i?j^^^(P) are independent of the direction of p. As we shall see in the next section, the 
attractive force around the Fermi surface picks the order parameter which is even with 
respect to the sign of for pairing. Therefore 

jrfi(p) = jm-p) (3-16) 

(3-17) 

jmn-jmn (3-18) 



and 

It follows from (3.7) that 
and 



In other words, the triggering source and the order parameter are symmetric with respect 
to simultaneous exchange of color and flavor indices. According to the antisymmetric and 
symmetric representations of the color group, we decompose further 

Brfi(p)=mi(n+xrfiin (3.20) 

with 
and 

xtm = xt%{p) = xttip) = xttip) 

with (f>'j]j^ (P) pertaining to the attractive channel and X%% (P) the repulsive channel. 
The color-flavor locked condensate for Nc = Nf corresponds, up to color-flavor-baryon 
number rotations, to 

Bf^iP) = B,{P)5%5% + B^{P)516%, (3.21) 
which gives rise to the minimum free energy of a homogeneous system below Tc- 
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4. The Di-Quark Propagator of QCD at High Baryon Density. 



The coefficient of T2{B) of (2.8) is the inverse of the di-quark propagator, which corre- 
sponds to the thermal diagrams in Fig. 2, when it is sandwiched by two order parameters. 
For total zero momentum and total zero Matsubara energy, the propagator is specified 
by the relative momentum, relative energy and color-flavor indices at both terminals. In 
writing, we have 

< P', s[,s'2; c[,c'2; f[,fi\V\P, s^, s^; ci, c^; /i, /s >= <^'^'^'^<^'^^'^<^/i///^/.<^P'pi^44,.is.(i') 

(4.1) 

where D represents the disconnected diagrams of the di-quark propagator, the flrst diagram 
of Fig. 2a, r is the sum of the proper vertices of the scattering between two quarks with 
zero total momentum and energy, and P = (p, —Vn) denotes an Euclidean four-momentum. 
The superscripts ci, C2, c'^, c'2 stand for color indices, the subscripts /i, /2, /(, /2 denote the 
flavor indices and the subscripts si, S2, s'l, s'2 and ti, ^2, ^'i, t'2 label the states above or below 
the Dirac sea. In order to facilitate the partial wave analysis, we found it convenient to 
associate the Dirac spinors u{'p) and v{p) to the vertex instead of the propagator. Thus 
the vertex in (4.1) is of the form 

r:'f4s::(^'i^) = ^{s'.^^ms'^. -^)r%fr\p'\p)Ua.{sum{s2, -p) (4.2) 

with the vertex function F^^q,^ given by conventional Feynman rules and U{s,p) = u{sp). 
The overall sign of (4.1) is flxed such that -D_|__|__(__|_ (P) h- > ^2_|_(p_^)2 as — > 0. 

The proper vertex function F satisfies the Schwinger-Dyson equation: 

+ m E (p'|p'')Ai.it...(P'')rS;r (^''1^) 

^ p//.ff ff f f^.f^'f ^ff 

(4.3) 

with F standing for two particle irreducible vertices. In terms of matrix notations, (4.1) 
and (4.3) take the form 

V^D- DTD (4.4) 

and 

F = f + TDF (4.5) 
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Combining (4.4) and (4.5), we find a simple relationship 

M = V-^ = D-^+T, (4.6) 

which generalizes the diagrammatics of a single particle propagator in terms of one particle 
irreducible self energy function. Upon decomposition of f into its symmetric (sextet color 
representation of SU (Nc)) and antisymmetric (anti-triplet color representation of SU {Nc)) 
components, 

I (4.7) 



2 

we obtain 

< P', S[, s'^\Ma\P, Si, S2 > +(5<^M^2C2 + 5c;c2^4ci) ^ p/^ ^/^^ s'^\Ms\P, Sl, S2 > 

(4.8) 

with 

< P\s',,s',\MAiS)\P,s,,s, >= 5p,pD-^,^^^^^^{P) + ^V'^!;!ls,^SP'\P)- (4-9) 

To the order of the approximation made in this work, the states below the Dirac sea do 
not contribute. We shall consider only the channels with all states above the Dirac sea 
and have the subscripts s's and t's suppressed. 

By calculating the matrix elements of M. between the order parameter (3.20), we 
obtain formally the quadratic term of the Ginzburg-Landau free energy functional with a 
single helicity 

= ^ E [< P'\-^a\P > Tr0t(P')0(P)+ < P'\Ms\P > TrxHP')xiP)i (4-10) 
^ P',p 

where we regard 4''f]j2 and x'f\% ^^^rix elements of NcNf x NcNf matrices (p and x- The 
symmetry factor | stems from the fact that we are only considering the condensate that 
consists of quarks of left helicity. With the HDL resummation, the irreducible vertex for 
the two quark scattering above the Dirac sea reads 

f ^(P'lP) ~ — - (l + —\ In — — ^-^ — - + higher partial wave components (4.11) 

12pp'\ Nc^ Wn-l^n'\ 



with 

2567r4V 2 ) n 
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and the function for two quarks propagating above Dirac sea is given by 

D~^{P) = [iUn- P + fJ'-'^iil^n,P)][-il^n - P + - ^(-it'n, p)] (4.12) 

with the single quark self-energy 

E(i.„,p).-^j!,^.„,„_g-^, (4.13) 

where mo — ^2^2^ denotes the Debye mass and is an infrared cutoff (|z/„| << Qc « jj)- 
As it was shown in [11], the weak coupling approximation in (4.11) and (4.13) is sufficient 
to determine the leading order of the pre-exponential factor of the transition temperature 
(1.1) or the gap energy. Since the color-symmetric channel does not carry the pairing 
modes, we may set 

<P'\Ms\P>^D-^{P)5p,p (4.14) 
with Dq{P) = the free di-quark thermal propagator. 

The spectrum of the operator Ma is given by the eigenvalue problem, 

J] < P\Ma\P' > u{P') = Eu{P) (4.15) 



I.e. 



Introducing 



we have 



D-\P)u{P) + ^Y. ^a{P\P')u{P') = Eu{P) (4.16) 

f{P)^-^Y.^A{P\P')u{P') (4.17) 

u{P) = ^ //^^ ^ (4.18) 

Substituting (4.18) back to (4.17), we find that the function f{P) satisfies a homogeneous 
Fredholm equation 

fiP) = ^Y.KE{P\P')f{P') (4.19) 

with the kernel given by 
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In order to demonstrate the equivalence with the methodology used in [10] , [11] to determine 
the transition temperature, we switch to the eigenvalue problem defined by the Fredholm 
equation 

KP) = ^Y.KE{P\P')h{P') (4.21) 

with the eigenvalues A^'s functions of the temperature, chemical potential and E. The eq. 
(4.19) implies that 

X^{T,n,E) = l (4.22) 

which determines E as function of T and For the pairing mode, we expect that £' = 
at T = Tc, i.e. A^(Tc,/i, 0) = 1, which is exactly the criterion for Tc in [11]. As T gets 
sufficiently close to Tc, we may write 

Ke{P\P') = K{P\P') + 5K{P\P') (4.23) 

with K{P\P') = -TAiP\P')D{P') and 5K{P\P') = -T a{P\P')D'^{P')E. Correspond- 
ingly, 

l = A-2(r,/i,0) + c(T, //)£;, (4.24) 

with the coefficient c(T, n) being calculated perturbatively. The eigenvalues of the kernel 
at £■ = and the transition temperature have been explored in [10], [11], and the one that 
satisfies the condition A^(Tc, 0) = 1 is given, in the vicinity of Tc, by 



+ ^) '"^ 7 b + 2^ - - + 4) i ^ 0(/)] , (4.25) 



A2 " 6^ 
where 

-mm. 



2567r3// 

and In ~ i, so that the 2nd and 3rd term inside the bracket are of the same order, 0{g). 
For |T - Tel « Tc, we have 

1 - 



^fl + ±Uln^ + ^ + ln2--(iV-l)(7r2 + 4)^^ 

A2(T,//,0) 37r4V nJV ec 16^ ' 



(4.27) 



37r^\ NJ ec Tc ' 

with ec = e|T=Tc- The constant pertaining the logarithm in the braket of (4.27), which 
contributes to the pre-exponential factor of the transition temperature (1.1) is subleading 
in the coefficient of T — Tc in comparison with In ^ and is dropped in the last step. 

For the rest of the calculation, we shall set T = Tc and maintain only the leading 
order term in g. Under such an approximation, the cigenfunction in (4.21) at = is 
given by [11] 

h(P) = —^=^^^sm( '^)e(6-\uJ) (4.28) 
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with 5 a cutoff energy such that ksTc « 5 « /i, and the eigenfunction u{P) in (4.18) 
that corresponds to the pairing mode reads 

u{P) = , ^'^-^^ (4.29) 
The choice of the constant in front of the sine of (4.28) normahzes u{P) according to 

By sandwiching (4.23) between h(P) and its adjoint, h(P) ~ ^^^(p'f]^)^ , and subsequently 
summing over P, we find that 

B C ec 

It follows from (4.24), (4.27) and (4.31) that 

7^42^c(T-Tc). (4.32) 

Note that, the summations over n in (4.30) and (4.31) are dominated by the terms with 
\n\ = 0{1) and the energy dependence of h{P) may be dropped to the leading order in g. 

Upon the decomposition of into its pairing mode and non-pairing modes, 

= 00 + 0', (4.33) 

with 

00 = \/6*w(P) (4.34) 

and 

-i^5^0j(P)0'(P) = O (4.35) 
the quadratic term of the Ginzburg-Landau functional (4.10) becomes 

r2(S) = 1^^4rc(T-Tc)Wt^+^{5^[^.^ + (p-/.)^]It[0't(P)0'(P)+x'(P)x(P)]}, 

(4.36) 

where is a constant NcNf x NcNf matrix, the function u{P) is given by (4.29) and 
we have made a zero-th order approximation (4.14), to the coefficients in front of the 
non-pairing modes. 
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5. The Quartic Terms. 



Because of the elimination of the pairing singularity in the diquark scattering channel 
embedded in the diagrams Q upon amputation, we may proceed with g — 0. The action 
becomes then quadratic and the path integral in the presence of the triggering term can 
be evaluated explicitly. 



Ignoring the Yang- Mills part and the fermionic states below the Dirac sea, the action 
(3.1) becomes 

SE{a, a) = Y.(-'^- +P- m,P^lp + Y^^hh iPr^)l-P^h,P + Jhh (^)«A,P«/L-p]' 
p p 

(5.1) 

where the repeated color-flavor indices are summed and the second sum, extends only 
half p-space. The triggering source for a given P — {p, —i^n)-, J'fl^i.P) element of a 

NcNf X NcNf matrix J(P), and the free energy shift is given by 



det 

W (J) = -kBTcJ2' In 

P det 



-il/n + P- 1^ J{P) 

J^{P) -iUn-p + IJ' 

-iUn + P- IJ' 

-il'n - P + P' 



(5.2) 



kBTc TrJt(P)J(P) ^ kBTc Tr[J^P) J (P)]^ ^ 



2 ^ J^^i + {P - P)'' 4 ^[ly^ + ip-fiy 



which corresponds to the sum of the ring diagrams in Fig. 3. Note the symmetry factor 
I of each term because of the single helicity pairing. Through a Legendre transformation, 
we obtain the quartic term of the Ginzburg- Landau free energy functional with a single 
helicity in terms of the order parameters 

T,{B) = ^ ^[^^ + {p- /^)TlV[St(P)S(P)]^ (5.3) 
p 

the coefficient of which is valid to the leading order of the coupling constant g. Upon the 
decomposition of the order parameter B into its color antisymmetric component and color 
symmetric one, eq.(3.20) and further decomposition of the color antisymmetric component 
into its pairing mode and non-pairing modes, (4.33), we arrive at 

UB) = ^J2^ul + ip- ^^ff^:r{[4iP)MP)? 
p 

+ 24{p)MP)[4iP)x{P) + xHp)Mp)] ^^'^^ 
+ 24{p)MP)[4iP)^'iP) + <i>'\P)UP)] + •••} 
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where • • • indicate terms of higher powers of x and cf)' . We would hke to emphasize here 
that the second and third terms of the right hand side of (5.4) arc in general nonzero. 
Indeed, for the color flavor locked condensate (3.21), at Nc — Nf — N 

i (5.5) 

we find that 

Tr(/>t(P)0(P)0t(p);^(P) = 

= ^^^'~^^^^~^^ \B,{P) - B^{P)\'[B^{P) - B^{P)Y[B,{P) + B^{P)]. ^^'^^ 

which is nonzero for N > 2. By combining (4.36) and (5.4), wc derive an expression for 
the Ginzburg-Landau free energy functional for a homogeneous system, 

r(B) = ^ Yl < p'\-Ma\p > TrcpUP')MP) 
P'P 

+ ^ Et^n + ip- ^ifm<p'Hp)<p'{p) + xHp)x{p)] 



2 

p 



+ ^ E(-n +ip- i^)rm4>Up)MP)f 



(5.7) 



p 

+ 2ct>Up)MP)[4>UP)xiP) + xHp)Mp)] 
+ 2(Pl{P)MP)[<PUP)<P'iP) + <P'Hp)MP)] + •••}, 

where a zeroth order approximation in g has been made for the coefficient in front of the 
non-pairing modes of the quadratic terms. Minimization with respect to x yields 

^h?.(p^ - --^n + (p - i^?mo{pWp)Mpmi + [MP)4^l{pmp)]7:f\ 
+ [up)4>l{p)up)]7:f\ + [up)4>l{p)upm\} 

and its contribution to the free energy is of the order of |0o|^- Below the transition 
temperature, 4>o{P) ~ (Tb — T)^ for Tc — T << Tc and the bulk condensate energy 
~ {Tc — T)^. We also notice the difference from the T = case, where the symmetric 
condensate is suppressed relatively to the antisymmetric one by a factor of order g [16]. 
On the other hand, the induced color-symmetric condensate, x ~ {Pc — T)2 and its 
contribution to the free energy ~ {Tc — T)^. The same is true for the non-pairing modes 
in the color-antisymmetric component, 4>'{P). Maintaining only the pairing mode, we may 
write 

0(P) ~ (f>o{P) = V6*w(P) (5.9) 
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where u{P) is the eigenfunction of the pairing mode, (4.29), and is a constant NcNf x 
NcNf matrix with the same symmetry of the indices as 0(-P). 

It follows from the expression (4.29) for u{P), (4.28) for h{P) and the eigenvalue (4.32) 

that 

r(;B) = n[-Tr*t* + ^Tr(*t*)2]. (5.10) 

with 

a=^4Tc(T-Tc) (5.11) 

and 

Including both helicities, we obtain the total Ginzburg-Landau free energy of a homoge- 
neous system with an even parity for the ground state 

TtotaiiBL^BR) = l][aWt^r + ^it(^t^)2j_ 

We conclude this section with some details for the realistic case of N^. — Nf — 3. With 
even parity the 9x9 matrix may be related to a 3 x 3 matrix via 

The trace of \E''s may be related to the trace of $'s. Explicitly, 

Tr*+* = 4tr$^$ 
T>(^rt^)2 ^ 2[(tr#$)^ + tr($^$)2]. 

Consequently, the Ginzburg-Landau free energy reduces to 

Vtotal = fi{4atr$t$ + 6[(tr$t$)2 ^ tr($t$)2]} ^^ -^^^ 

where a and b are given in (5.11) and (5.12) at weak coupling with a < for T < Tc and 
b > 0. The minimum free energy below Tc corresponds to 

$t$=M, (5.17) 

i.e. $ is proportional to an unitary matrix. By a flavor and baryon number rotation, it 
can always be brought into an unit matrix, which leads to the standard expression of the 
color-flavor locked condensate with 
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The corresponding condensate energy is 



{TtotalUn = -6-n = -^/c^T^(^^^ j n. (5.19) 

We have thus proven that the color-flavor locked condensate remains energetically favored 
near Tg. Using (5.8) and (5.9), the induced sextet order parameter reads 

and it contributes to the condensate energy a term 

_ 558C(5) .j.2t2( Tc-T Y 

which is smaller than (5.19) by a factor -2:^^^. The existence of the sextet order param- 
eter is evident from the composition rule 3x3x3 = 3-|-3-|-6-|-15of the color SU (3) 
representations and its conjugate for the flavor SU{3) representations. The sextet correc- 
tion (5.21) may be compared with the correction due to the constraint of color neutrality 
discussed in [15]. The latter is suppressed by a factor of (^^p-)^ In relative to (5.19). 

If > (^^^)^ln the sextet correction will be larger in magnitude. 



6. Kinetic Energy and Gauge Coupling. 



An important application of the Ginzburg-Landau theory lies in systems which include 
inhomogeneous condensates, such as domain walls and vortex fllaments. To incorporate 
the spatial variation of the condensate, a nonzero total momentum is introduced to the 
di-quark propagator. As we shall see below, the characteristic momentum corresponds to 
the inverse of the coherence length, which is much longer than the thermal wavelength 
of the system and diverges at T^. Therefore, the momentum dependence of the di-quark 
propagator will be treated perturbatively while the dependence of the quartic coeflicient 
will enter only through the constraint of the overall momentum conservation. 

A non-zero total momentum k will modify both the disconnected and the connected 
part of the di-quark propagator. The spectral problem (4.16) is replaced by 
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where we have appended the subscript k to indicate the dependence on the total momentum 
while the relative 4-Euclidean momentum remains denoted by P. Thus the 4-Euclidean 

momentum of each quark in a pair is (| ±p, =F'^n)- The homogeneous case considered 

in section 3 corresponds to the special case k = 0. As it was shown in section 3, the 
quark self-energy may be omitted to the leading order of the interaction, except in the 
determination of the transition temperature, DZ^{P) ^ {ivn — p+ + ii){—Wn — P- + /^) 

with p-i- = || ±p|. The Fredholm equation (4.21) becomes 



with 



The expansion (4.23) is supplemented by more terms 

^E,k(P\P') = K{P\P') + 5K{P\P') + 5'K{P\P') + 5"K{P\P') (6.4) 

with 



{P - k- - {kpf 



(6.5) 



4p ul+{p-ii)^l 



while the term S" K{P\P') arises entirely due to the momentum correction of the vertex 
function. The corresponding eigenvalue condition (4.24) is replaced by 

l = X-^{T,i^,0) + c{T,iJ,)E + c'{T,iJ,)k^ (6.6) 

where c'{T,n) is computed perturbatively, like c{T, /j,) of (4.31). The first order perturba- 
tion yields the result 

^ -24.^J|?iln^ = -r^^^' 

This contribution arises exclusively from the second term inside the bracket of (6.5), 
while the third term contributes a term smaller by a factor whose order is less than 
(^tP')^ fciT"- '^^^ second order perturbation of the first term inside the bracket of 
(6.5) contributes to c'(T, //) a term that is smaller than (6.7) by a factor of order g. Fi- 
nally lets comment that the contribution of the term 5" K{P\P') is analyzed in Appendix 
B, where we show that its contribution to c'{Tc,iJ,) is comparable with that of the third 
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term inside the bracket of (6.5) and therefore can be ignored. Finally, the eigenvalue of 
the pairing mode with a small total momentum becomes 

E^^4Tc(T-Tc) + ^P. (6.8) 



Therefore, the typical momentum of the condensate is ~ ksTcy whose inverse 

corresponds to the coherence length. For a condensate consisting of superpositions of 
different ^'s, i.e. 

B{P) = V6^^U^j;{P) (6.9) 

k 

the Ginzburg-Landau functional becomes 

r = ^ y (fc2 + a)Tr*t*r + ^6 V mI *t ^.-^r r '^r ' (6.10) 

2 ^ fe 4 kl k2 fcl+fc2,fcl +k2 ^ ' 

k k\k2k\ k2 

By transforming the condensate from the momentum space into the coordinate space via 

*(-"0 = ;^E*.-^^''^ (6-11) 
k 

and including both helicities with an even parity ground state we obtain the conventional 
form of the Ginzburg-Landau free energy 

Ttotal= J dVjltV^tv* + aMt$ + i5Tr(*t^)2 . (5.12) 

It is instuctive to compare the Ginzburg-Landau free energy (6.12) with that of the ordinary 
BCS model. On setting vp = 1 for the latter and scaling 1— ^^/2kF'^j^, the free energy 
as a polynomial in \E'^ contains identical coefficients as (3.21). 



It remains to couple the order parameter to the classical Yang-Mills color gauge field 
and to the electromagnetic field. While the coupling structure is dictated by diagrams, we 
shall follow a different route by replacing the gradient in (6.12) with the gauge covariant 
derivative. As the final result ought to be gauge invariant, such an approach is unambigious. 
Let A — A^T^ denote the classical vector potential of the SU (Nc) color gauge field and A 
that of the electromagnetic field, of which the associated U{1) rotation is a subgroup of 
the fiavor SU{Nf) group. The gauge covariant derivative acting on a single quark field is 

{Di;'^) f = Vrf- igA^^'i^f - ieq^ri (6-13) 
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with qf being the charge number of the /-th flavor. Prom here on, only the repeated color 
indices in (6.13) are summed. Since the order parameter behaves like a direct product of 
two quark fields, the gauge co variant derivative of ^ reads 

[5^)11 = v^ii - 19!'^^'^'^%% - igA^^'^tn - '<^h + ^if.Un.t (6-14) 

By replacing V^' of (6.12) with and including the energy of the pure Yang-Mills part, 
we obtain the gauge invariant Ginzburg-Landau free energy functional of an inhomogeneous 
superconductor 



total 



/r 1 1 1 "1 

d^f -F^jF^j + -{^xAf + Tr(^^')t(I)^r) + aWt$ + -6T>(^'t^)2 ^g^^^^ 
. Z Z Z - 



It is easy to verify that the free energy functional is invariant under the local gauge trans- 
formations 

/1 72 71 J2 

Ai-^uAu^-\ — uWu^ (6.16) 

- 1 - 
A Va 

e 

The equation of motion follows from the free energy functional by applying the variational 
principle. 

Finally, let us explore some details for the realistic case of = Nf = 3. It follows 
from (6.13) that 

m^'^-)f,f, = e'^^^-^ef^f^fiD^)} (6.17) 

with 

{p^)) = ^^)-igA ^-ieQfA^ (6.18) 

where A = tI't' with T' = — t'* being the generator of the 3 representation, and Qf = 
?/i + 9/2 //1/2 represent a cyclic permutation of 1, 2, 3. Arranging the flavor index in 
the conventional order of u, d, s we find the diagonal electric charge matrix 

Q^dtag{^,~,-^) = ^T\ (6.19) 

Here we have adapted an expression of that differs from the standard one by a cyclic 
permutation of rows (columns), in which = -^^diag{2, —1, —1). The equation (6.18) 
takes the matrix form ^ 

D^ = V^- igA^ - leqA^f^ (6.20) 

with q = —-^e. On noting that $ = ^>o + T'^^, we have the decomposition 

= {D^)o + f^D^i (6.21) 
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with (-D$)o = |tr(^^>) and D^i = 2trT'D$, and its contribution to the free energy (6.15) 

tr(^$)t(^$) = 3{D^)*{D^)o + ^{D^)1{D^)i (6.22) 

Now we divide the field variables into two groups, the first group contains $o, ^8, ^8,^, 
while the second one includes A'' with / ^ 0, 8. A closer inspection of the structure of 
L>$ reveals the absence terms in the free energy (6.15) that are linear in the variables of 
the second group. Therefore a solution of the equation of motion exists in which only the 
fields of the first group acquire non-zero values. In what follows, we shall focus on this 
special case. 

By setting = A'^ = for I ^ 0,8 and transforming the remaining variables 

= $o + ^$^ x=^(*o-^$^) 

' 2V3 (6.23) 

F = videos 6* + ^sin 6* V = -i'^sin6' + ^cos6', 

with tan 6* = we arrive at the following expression for the Ginzburg-Landau free 

energy functional (6.15) 

Ttotai = [ c^V(J(V X Vf + J(V xVf + \{V- 2iAV)0|2 + |(V + iXV)xf 

J ^ (6.24) 

+ 4a(|</>p + IxP) + + \x\^f + |</>r + \\x\% 



with A = i^j^s^ 9^ + The transformation of the gauge potentials in (6.23) is nothing 
but the transformation performed in [3], [21], [22] that selects out the unbroken U{V) 
gauge potential V in the color-flavor locked phase. The V held remains decoupled from 
the condensate in the presence of a nonzero and can be set to zero as a special solution. 
The condition for the color-flavor locking (that the 3x3 matrix is proportional to an 
unitary matrix), 

|X|'=2|0|2 (6.25) 

which corresponds to the free energy minimum in the bulk, cannot be imposed everywhere 
because of the nonzero values of V0, Vx and V , as it is evident from the example of 
a vortex filament and a domain wall which are analyzed in the appendix C. In another 
words, in the presence of an inhomogeneity and of a non zero gauge potential, the unlocked 
condensate, which is an octet, under a simultaneous color-flavor rotation mixes with the 
locked condensate. This situation may be compared with the s — d wave mixture in the 
core region of a vortex fllament of a cuprite superconductor [23] . 
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7. Concluding Remarks. 



In this work, wc have derived the Ginzburg-Landau free energy functional of the color 
superconductivity from the fundamental QCD action by means of normal phase thermal 
diagrams. We have also addressed the issue of the induced color symmetric component 
of the condensate near Tc and the interplay between the color-flavor locked condensate 
and the unlocked one in the presence of a gauge potential. Though the expressions of the 
Ginzburg-Landau coefficients were derived in the weak coupling approximation, which is 
valid in the ultra-high baryon density regime, we expect that many qualitative features, 
such as the color-flavor locking, the mixture of the color-symmetric condensate and the 
coupling among the different types of condensates and the gauge potential, to survive in 
the regime of realistic baryon density, e.g., the baryon density achieved at RHIC or that 
inside of a neutron star. 

As the chemical potential is lowered, higher order corrections become more important. 
There are two types of higher order terms, which constitute a hierarchy of smallness, 
those which result from higher powers of the coupling constant g and those that arise 
from the curvature of the density of states at the Fermi level. It follows from (1) that 
g ~ \n.~^ /i/kBTc, while the latter corrections come in powers of {IibTc / ■ Therefore 
the corrections due to the higher order thermal diagrams are far more important. In view 
of the 27r suppression of the integral over loop momenta inside higher order diagrams 
(before the combinatorial factor dominating the magnitude), the numerical values of the 
Ginzburg-Landau coeflBcients we derived may still serve as a reasonable approximation for 
a realistic high baryon density. 
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In this appendix, we shall illustrate our methodology by rederiving the Ginzburg- 
Landau free energy for the BCS model of the system of nonrelativistic electrons. The 
model Hamiltonian is given by 

H = Y1 ^P^ls^P^ ^^^^ ~ \^p\Wi<^D - \^^\)ai^al^^a_^iap'^, (9.1) 

ps p,p' 

2 

where = 2^ ~ A*? ^ the electron mass and fj, the chemical potential. The subscript s 
labels the spin state, ^ > is the pairing coupling and cud denotes the Debye frequency. 
The path integral expression of the free energy at a temperature T reads 

'^''p ( " fcS") " ^/ n^^p^^^p^^"^"^"''^ (9-2) 

ps 

where the Euclidean action 

Ssia, a) = '^{-ivn + ^p)apsaps 

Ps 

~ Wl ^ ^^^^ ~ \ip\)^{^D - \ip'\)o^p^^a^p_iO'P'_iap!^^5n++n-,n'^+n'_-, 
P+P-,P'_^,P'_ 

(9.3) 

where P = {p, — fn) is the Euclidean four momentum and P± — (±p, ^v). The Matsubara 

energy = 2'KnkBT with n = ±i,±|, Following the steps outlined in section 2, we 

add to Se a triggering term 

= Y.^{uJD - M)[r{P)a.p^ap^ + J(P)at^aLpJ, (9.4) 
p 

which gives rise to an increment W{J) of the free energy and induces an order parameter 

The Legendre transformation oiW{J), 

V{B) = T2 + T^ + ... (9.6) 

is the Ginzburg-Landau functional and determines the order parameter at equilibrium 
when the triggering term (9.4) is removed, with 

T2^^B*{P) <P\M\P' > B{P') (9.7) 
p,P' 
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and 

r4 = ^ B*{Pi)B*{P2)<PiP2\Q\P[P2>B{P[)B{P!,). (9.8) 

The coefficient < P\M\P' > is the inverse di-electron propagator, given by 

< P\M\P' >= {vl + il)5pp, - j^e{u^D - \ipmu^D - kp'l), (9.9) 

where the second term on the right hand side represents the leading order 2PI vertex in 
accordance with the Hamiltonian (9.1). 



The eigenvalues and eigenfunctions of M. are determined by 

9 



i^n + iiMP) - iLoiu^D - mJ2^{u^^ - i^p'MP') = (9-10) 



On writing 

A = 

pn 



P' 

we have 

Substituting (9.12) back to (9.11), we obtain the secular equation for E 

The pairing mode corresponds to the eigenvalue that vanishes at Tc- For T, sufficiently 
close to Tc, we may expand the r. h. s. of (9.13) to first order in E, i.e. 



(9.14) 



gKp / 2UD \ ((,{6)Kf 



with kp the Fermi momentum, vp the Fermi velocity and 7 the Euler constant. At the 
transition temperature T = Tc the following relation holds, 

l=^flnzl^+T) (9.15) 



2n'^vp \ TrkpTc 
and E = according to (9.14). In the neighbourhood of Tc, we have 



E=^klTc{T-Tc). (9.16) 
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All other modes are nonpairing with eigenvalues Ei = i^n'^^p + ^id)- Expanding the order 
parameter according to the eigenfunctions of A4, i.e. 

B{P) ^-^uiP) + B'{P) (9.17) 

with u{P) the eigenfunction of the pairing mode and ^ J^^p tt'''(P)S'(P) = 0, we find 

T2 = E^*^ + ^(i/2 + ^/)B'%P)B'{P). (9.18) 
p 

To calculate the quartic term, we switch ofi^ the interaction and obtain the quartic term of 
J in W{J). Then we amputate it with four di-electron propagators. The result reads 



r4 = 



1*1 

2pU 



(9.19) 



where the terms in ... contain nonpairing modes and are at most ~ \'^\^ as T 1— > Tc- 
Neglecting these modes, we find 



r(*) = Q 



7C(3) 



■(T-Tc)|*r + 



7C(3)4 



(9.20) 



— * 

If a small momentum k is introduced, only the eigenvalue of the pairing mode is 
modified by an amount comparable to itself. In order to determine E in this case, we 
observe that (9.13) is replaced by 



^ {iUn - ^p+){-Wn -^p_)-E' 



(9.21) 



where p± = ^ ±p. For the pairing mode, we expand the r. h. s. to the first order in E 

— * 

and to the second order in k. After some manipulations, we end up with 



7C(3) Tc 6 ^ 



(9.22) 



On writing \E'^ug(P) for the order parameter at a nozero k with Uj^{P) being the pairing 

— * 

mode at k ^ 0, the free energy functional becomes 



^ ^ -6 ^ 7C(3) Tc 



1^ 2 , S7TW{kBTc)\ |4 

^ 7C(3)4 ' ^' . 



(9.23) 
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By summing over k and introducing 




3 J (27r)3 ^ 

we obtain, for an inhomogeneous order parameter, that 



-i 



1. 



2m 2 



with coefficients 



and 



a = 



-klTc{T-Tc) 



367rV {ksTcf 



7C(3) //^^ ' 
which is in agreement with the result obtained by Gorkov [13]. 



(9.24) 



(9.25) 

(9.26) 
(9.27) 
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In order to determine the vertex contribution to the coefficient c'(Tc, //) in eq. (6.6), 
we approximate ^[P\P') by considering only the one magnetic gluon exchange, with 
both incoming and outgoing quark momenta near the Fermi surface, i. e. 

^Ark(P\P') - -^(1 + ^)i^M(g,a;)^(pV)7»«(P+)^(p'-)7i«(p-)(% - mj) (10-1) 
where u{p) is a solution of the Dirac equation (74^ — i-f • p)u{p) = u^{p)u{p) = 1, and 

P± = ^±P, p'± = |±p', q = p-^ (10.2) 

with p p' and Dm{(1, uj) contains the HDL resummation. With the aid of the Dirac 
equation, (10.1) can be reduced to 



^A,mP') - -^(1 + ^)DM{g,u;)[V - ^l±.J±}^^^s] (10.3) 



where 

V = - v) {^^)au{p+)u\^_)au{p-) 
S =u\'^^)u{p^)u\^_)u{p_) 
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with a = and r the 2x2 Pauh matrices. Since u{p) is also an eigenstate of 

15= (^^ J we have 

with Tp(f){p) = (t){p). Using the identity {T^)a/3{r^)pX = '^^aX^Pp - ^ap^pX we find that 

V = S- 2S', S' = c^\^^)c^{p-)4.\^_)<^{p+). (10.6) 
In the collinear limit, g i— > 0,p± ^ p'^-, we have S ^ 1 and 

S' ^ ^(1 +P+P-) ~ ^[1 - (p^)^]. (10.7) 

Therefore, a net momentum k shifts the kernel K{P\P') by an amount 

5"K{P\P') ^ -K{P\P')^[l - ip'kf]. (10.8) 

Comparing 5" K{P\P') with the second term of 5' K{P\P') in (6.5), we notice that although 
the former differs by a factor of + (p — /^)^ in the denominator in comparison with the 
latter, it compensates by the factor of p^ ~ //^ which appears in the denominator. As the 
sum over Matsubara energies will pick up the infrared logarithm. In |, we expect that the 
contributuion of 5"K{P\P') to c'{Tc,iJ,) is suppressed relative to that of of 5'K{P\P') by 
a factor of order 

(^)'ln^«l (10.9) 

\ 11 J kBJ-c 

and therefore can be ignored. 
11. Appendix C. 



By applying the variational principle to the free energy functional we obtain the 
equation of motion for the special case oi — Ai — for / 0, 8 and V = 0, 

V X (V X F) = J 

- (V - 2iXVf(f) + 4a(f) + 26(2|(/.p + |x|^)0 = (nj) 

- (V + iXVfx + 4ax + 2b{\cP\^ + ^ |xP)x = 
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with the current vector given by 



J = iX 



20*(V-2iAy)0 + X*(V + iAy)x +•••, (11.2) 



and the ■ ■ ■ indicate the complex conjugate. For a vortex filament parallel to ^, we make 
the following ansatz for the solution in cylindrical coordinates 

= ^p)e2~^ X = ^(p)e-~^ V = V{p)e (11.3) 

The set of equations (11-1) is reduced to a set of ordinary non-linear differential equations 

dp p dp p 

- ^ - - ^ + 4(- - XV fu + 4au + 2b{2u' + v^)u = (11.4) 
dp^ p dp p 

-— r -+ XVfv + Aav + 2h(u^ + -v^)v = 

ap^ pdp p 2 

If the color and the flavor were locked everywhere, we would have v — \/2u. The last two 
equations would then imply that 

F=f (11.5) 

Ap 

which is a pure gauge for p 7^ 0, as the total flux has been flxed by the phase of (j) and x 
in (11.3). A complete color-flavor locking would squeeze all flux in the symmetry axis of 
the filament and would make the magnetic field energy per length infinite. 

The next example is the equilibrium between the super phase and the normal phase 
under an external V-field, the classical problem of the surface tension of the interface which 
was worked out by Ginzburg and Landau for electronic superconductors [24]. With the 
interface parallel to the y — Q plane, the solution ansatz reads 4> — u{x), x — 'f^(^) cind 
V = V{x)(, where u,v are real and positive. The set of equations (11.1) reduces to 



dx'^ 
d^u 



+ 4:X'^V'^u + 4au + 2b{2u'^ + v^)u = (11-6) 



dx"^ 



+ AV^-y + 4.av + 2h{u^ + -v'^)v = 



dx^ ' 2 

An inspection of the last two equations shows that a complete color-flavor locking, v = V2u 
can not implemented unless V — 0. 
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(a) (b) 

Figure 1 

The continuation from the normal phase to the super phase 

a) The continuation path on T — J plane 

b) The order parameter versus temperature 




(a) (b) 

Figure 2 

a) Diagram of a quark pair propagator. The thick black line indicates the full quark 
pair propagator. 
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b) Dyson-Schwinger equation for the proper vertex of diquark scattering, where the 
vertically shaded vertex stands for two quark irreducible diagrams. 




Figure 3 

The diagrammatic expansion of W(J) at g=0. 
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